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Abstract. 

If we assume that we live in the center of a spherical inhomogeneous universe, we can 
explain the apparent accelerating expansion of the universe without introducing the unknown 
dark energy or modifying gravitational theory. Direct measurement of the cosmic acceleration 
can be a powerful tool in distinguishing ACDM and the inhomogeneous models. If ACDM 
is the correct model, we have shown that DECIGO/BBO has sufficient ability to detect the 
positive redshift drift of the source by observing gravitational waves from neutron star binaries 
for 5-10 years. This enables us to rule out any Lemaitre-Tolman-Bondi (LTB) void model with 
monotonically increasing density profile. Furthermore, by detecting the positive redshift drift 
at z ~ 0, we can even rule out generic LTB models unless we allow unrealistically steep density 
gradient at z ~ 0. We also show that the measurement accuracy is slightly improved when we 
consider the joint search of DECIGO/BBO and the third generation Einstein Telescope. This 
test can be performed with GW observations alone without any reference to electromagnetic 
observations. 



1. Introduction 

When we assume that our universe is homogeneous and isotropic, current cosmological 
observations (e.g. type la supernovae (SNe) [1]) indicate that the cosmic expansion is 
accelerating. Once we allow the possibility that our universe has cosmological scale spherical 
inhomogeneity with the observer at the center, the observations can be explained without 
introducing the unknown dark energy or alternative theories of gravity. In such models, 
the Copernican Principle is apparently violated and the cosmic expansion is not necessarily 
accelerating. Therefore, future direct detection of the cosmic acceleration is very useful in 
distinguishing them. The direct detection of the acceleration of the universe provides not only 
a key to solve the dark energy problem, but also a critical test of the Copernican Principle. 

The simplest example of the inhomogeneous model is the Lemaitre-Tolman-Bondi (LTB) 
spacetime which is a spherically symmetric dust solution of the Einstein Equations. The metric 
is given as 




(1) 



where R(t, r) and k(r) are arbitrary functions. These functions are related to the density of 
the dust p(r) via Einstein Equations. If we live at the center of the LTB spacetime with 
a Gpc-scale void, the apparent cosmic acceleration can be explained. The LTB model has 
been partially tested with the cosmological observations like, the cosmic microwave background, 
baryon acoustic oscillations, the kinetic Sunyaev-Zeldovich effect etc., [2-14], but it has not 
been completely ruled out yet. For example, Ref. [13] claimed that LTB void models are in 
conflict with current observations with nearly scale invariant primordial spectrum, but Nadathur 
and Sarkar [14] claim that that observational results can be explained by assuming different 
primordial spectrum ansatz. Therefore, the observations which are unaffected by the primordial 
information are very crucial to test a wide class of LTB models. One of such observations is the 
redshift-distance relation of type-la SNe. However, it has been shown that one can construct the 
LTB void model that exactly reproduces the redshift-distance relation in ACDM [15]. Therefore, 
we need other observations that do not depend on the primordial information. 

Redshift-drift measurement is the one that meets our demands [16] (see also Refs. [15, 17, 18]). 
Redshift drift is the time evolution of the redshift due to the cosmic acceleration, hence its 
detection means the direct measurement of the acceleration of the cosmic expansion. 

In the Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime, the redshift drift is given 

as 



where At Q denotes the observation period, and Hq and H(z) are the Hubble parameter at 
present and at redshift z [19], respectively. In the ACDM universe, Atz is positive in the range 
z = — 2 [17]. On the other hand, Atz in LTB spacetime obeys the following differential 
equation [18]; 



where R(t, r) is an arbitrary function. Recently, Yoo et al. [18] have shown that when the matter 
density is monotonically increasing, the right hand side of the above equation is negative. By 
combining this with Afz| 2= o = 0, we can show that Atz must be negative at any z in this model. 
Furthermore, they have shown that for any LTB density profile, dAtz/dz < for z <C 1 unless 
we allow unrealistically steep density profile at z ~ 0. Therefore, it is crucial to measure the 
sign of At z at z < 2 in order to distinguish the ACDM and LTB models. 

The order of magnitude of Atz is roughly given as the cosmic age divided by the observation 
time, hence A^z ~ 10~ 10 for lyr observation. It is this tiny value that makes it difficult to 
measure Atz with current technology. Recently, Quartin and Amendola [17] have shown that 
by measuring the shift of the Lyman a forest of quasar spectrum at z = 2 — 5 with the proposed 
E-ELT instrument CODEX [20] for 10 yrs, it will be possible to distinguish ACDM and typical 
LTB void models. However, CODEX would not be able to measure Atz at low z since Lya 
forest can be measured from ground only at z > 1.7 [20]. Hence, they can only test typical LTB 
models but not generic ones. 

In this paper, we estimate how accurately we can measure the redshift drift with future 
gravitational wave (GW) interferometers. It seems that DECIGO [21,22] and BBO [23] are the 
only proposed detectors that can measure Atz at z < 2. We consider neutron star (NS) binaries 
as GW sources, which are often called as the standard sirens and can be unique tools to probe 
the cosmic expansion [24-26]. When the expansion is accelerating, we may find an additional 
phase shift in gravitational waveforms [21]. We assume that ACDM is the correct model and 
estimate whether we can tell the positivity of the redshift drift at low z with GW observations. 

Throughout this paper, we use the unit G = c = 1. 
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2. Correction in gravitational waveform due to the redshift drift 

Let us first derive the correction in the gravitational waveform phase due to the accelerating 
expansion of the universe. We here consider a binary consisting of two bodies with masses mi 
and 7TT-2- We define the time to coalescence measured in the observer frame as At = t c — t with 
t c representing the coalescence time. This At includes the effect of cosmic acceleration. On the 
other hand, we denote the time to coalescence measured in the source frame as At e and define 
AT = (1 + z c )At e where z c is the source redshift at coalescence. Then, the relation between At 
and AT is [21,27] 

At = AT + X(z c )AT 2 , (4) 
where X(z) is the acceleration parameter defined as 

X(z) = -° 1- }> . (5) 



2 V (l + z)H . 
Notice that X(z) is related to the redshift drift AfZ as 

A t z = 2(l + z)At X{z). (6) 



By using the stationary phase approximation [28], the waveform in the Fourier domain can be 
expressed as 

M/) = e^ c{/) M/)| noaccel , (7) 

where 

*aoc(/) = -2irfX(z c )AT(f) 2 (8) 

and /*(/)! 

no accel corresponds to the gravitational waveform in the Fourier domain without cosmic 
acceleration explained in the Appendix. The leading AT(f) is given as [28] 



AT(f) = 5(8ttMJ)- 8 / 3 M z , (9) 

where A4 Z = M(l + z c )i] 3 ^ denotes the redshifted chirp mass with M = (mi + m-i) and 
rj = mim^lM 2 representing the total mass and the symmetric mass ratio, respectively. With 
this at hand, ^ / a cc(/) is given as [21,27] 



= X(zc)M z x~ A , (10) 

where x = (TrA4 z f) 2 ^ 3 . A term proportional to *$>N(f)x n represents the ro-th post-Newtonian 
(PN) order relative to the leading ^jv(/) = ifg (^-Mzf)' 5 ^ 3 , hence this is "-4PN" correction. 
We used the restricted-2PN waveform including spin-orbit coupling at 1.5PN order , where 
"restricted" means that we only take the leading Newtonian quadrupole contribution for the 
amplitude and neglect the ones from higher harmonics. 



3. Numerical Setups 

The measurement accuracy on Atz has been estimated mainly using ultimate DECIGO (which is 
three orders of magnitude more sensitive than DECIGO) [21, 27]. Especially, Ref. [27] used this 
measurement accuracy on the acceleration parameter X(z) to estimate the ones on cosmological 
parameters. We improve their analyses in the following way and apply the result to probe 
the inhomogeneity of the universe: (I) We use DECIGO/BBO (and not ultimate DECIGO), 
including the confusion noises from white dwarf (WD) binaries. (II) Rather than sky-averaged 
analysis, we perform Monte Carlo simulations by randomly distributing the directions and 
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Figure 1. The root noise spectral density of DECIGO/BBO with (solid) and without (dashed) 
WD confusion noises, and the one of ET (dotted). 



orientations of sources. (Ill) We include the spin-orbit coupling into binary parameters. (IV) 
We use the merger rate that reflects the star formation history. 
We take the binary parameters as 



9 i = (lnM z ,lnri,P,t c ,(l>c,D L ,9 s ,(j) S ,6 L ,(t) L ,X H ) 



(11) 



Here, (3, (fi c and Dl represent the spin-orbit coupling parameter, the coalescence phase and the 
luminosity distance, respectively. (0s,4>s) are the direction of the source in barycentric frame 
which is tied to the ecliptic and centered in the solar system barycenter (see Fig. 1 of Ref . [29] ) , 
and (#Lj 4>l) are the orientation of the source orbital axis in the same frame. We have introduced 
a new parameter Xh which is defined as Xh = X(z)/Hq. Xh is related to AfZ through Eq. (6). 

We estimate how accurately we can measure binary parameters 9 l using Fisher analysis. 
Assuming that the detector noise is stationary and Gaussian, the measurement accuracy is 

/-_ \ l/2 

given as A.9' 1 = IT J [28], where, Tij is defined as 
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(12) 



Here, 50 % = 6 % — 0\ rue with #tme representing the true values of binary parameters and (9) is 
the prior information. Tij is the Fisher matrix which is defined as 



4Re 



/fl " , f h*(f)hf) 
/ in ; S n (f) 



(13) 



where S n (f) denotes the noise spectrum of the detector. In this paper, we use the instrumental 
noise spectrum of BBO when we perform Fisher analyses and assume that DECIGO also has 
the same sensitivity as BBO. The noise curves are shown in Fig. 1. The actual expression of the 
total noise spectrum is given in Eq. (36) of Ref. [30] . 
/in and /fi n in Eq. (13) are given as [30] 



/ in = (256/5)- 3 / 8 7T- 1 A^; 5 / 8 At- 3 / 8 , / fin = 100Hz 



(14) 



respectively. f m represents the frequency at At before coalescence and /fi n denotes the higher 
cutoff frequency of DECIGO/BBO. For fiducial values, we set m\ = rri2 = IAMq and take 
tc = 4>c = P = 0- Since the dimensionless spin parameter is less than 1, we adopt the prior as 
\/3\ < 9.4 [31]. We assume that a flat ACDM is the correct model. In the following computations, 
we set the fiducial values as Hq = 70km/s/Mpc and the cosmological parameters as O m = 0.3 
and Oa = 0.7. 

The number of binaries AN(z) that exists in each bin with size 5z = 0.1 is estimated as [32] 

AN(z) = 4vr [a r(z)] 2 h(z)(dr/dz)6zAt , (15) 

where 

a r(z) = [ dz'/H(z'), dr/dz = {(1 + z)H(z)}- 1 (16) 
Jo 

with the Hubble parameter at redshift z given as 



H( z ) = H ^n m (i + z)3 + n A . (17) 

Here, ao, f(z), and r each represents current scale factor, comoving distance and proper look 
back time, respectively. h(z) = hoR(z) shows the NS/NS merger rate per unit comoving volume 
per unit proper time, where we assume the merger rate today as no = 10 -6 Mpc~ 3 yr -1 [33] 
and the merger rate evolution against z as 

( 1 + 2z (z < 1) 
R(z) = { |(5 -z) (1<2<5) (18) 
I (z > 5), 

which is based on the current observation of star formation history [34]. 

Following Refs. [29-31], we randomly generate 10 4 sets of (9s, 4>Si 4>l) for each fiducial 

/-_ ,N 1/2 

redshift Zf. = 0.1k + 0.05 (k = 0, 1, • • •). Then, we calculate IT j for each set and take the 



average to yield 



1/2 



. The measurement accuracy at each z^ is estimated as 



A9 l = N-~t AN(z k 



f- 1 



1/2 



(19) 



where A^; n t = 8 shows the number of effective interferometers. DECIGO/BBO consists of four 
clusters of triangular detectors (see e.g. Fig. 2 of Ref. [30] for the proposed configurations of 
DECIGO/BBO) and for simplicity, we assume that all the clusters are placed on the same site. 
We also use the sky-averaged analysis to see how the new effects that we considered in this paper 
affect our results. 



4. Results 

In Fig. 2, we show the measurement accuracies of A^z using DECIGO/BBO for 5 yr observations. 
We see that the positivity of A t z can be marginally detected at redshift around z ~ 0.5. Figure 3 
is same as Fig. 2 but for 10 yr observations. If this observational period is realized, AfZ > can 
be detected with 3-cr confidence level. 

Next, we show the effects of including f3 and WD confusion noises. Figure 4 shows 
A(Afz)/ AfZ for the sky-averaged analyses with 5 yr observations where the sky-averaged 
waveform is given in Eq. (A. 12) 1 We see that the spin-orbit coupling (3 affects the result 



We have assumed multiple detectors, but in principle, sky-averaged analysis can be performed with a single 
detector. 




Figure 2. The measurement accuracies of AfZ using DECIGO/BBO for 5 yr observations. 




Figure 3. The same as Fig. 2 but the observation period changed to 10 years. 



on lower z side. This is because the redshifted mass M z is lower for lower z source, which leads 
to smaller orbital velocity v at a given frequency / and larger degeneracies between ft and other 
binary parameters. On the other hand, WD confusion noise affect higher z side. This is because 
higher z source leads to GW signal with lower frequency where the effect of the confusion noise 
is larger. 

Similarly, let us compare the results using sky-averaged analysis and Monte Carlo simulation. 
Figure 5 shows the former in solid and the latter in dashed lines. Also, sky-averaged analysis 
without taking ft nor WD confusion is shown in dotted line. By comparing the dashed curve 
with this dotted one, we see that taking ft, WD confusion noise, directions and orientations of 
sources into account deteriorate the determination accuracies of At z by a factor of a few. Since 
5 yr observation can only marginally detect the positivity of AfZ, this difference is crucial in 
detecting A t z > 0. 
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Figure 4. A(Atz) / AfZ using sky-averaged analysis for 5 yr observations, taking both j3 and 
WD confusion noise into account (solid), without (3 (dashed) or without WD confusion noise 
(dotted). 
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Figure 5. Solid curve is same as the one in Fig. 4 while dotted curve does not include (5 nor 
WD confusion noise. Dashed one shows the result with the Monte Carlo simulation. 

5. Conclusions and Discussions 

In this paper, we have estimated how accurately we can directly measure the cosmic acceleration 
using future space-borne GW interferometers such as DECIGO and BBO. If we assume that 
the ACDM model is correct, we have shown that we will be able to measure the positivity of 
the redshift drift with 5-10 yr observations, which enables us to rule out any LTB void model 
with monotonically increasing density profile. Furthermore, since 10 yr observation allows us to 
detect the positive redshift drift at z ~ 0.05, it seems that we can even rule out generic LTB void 
models unless we allow unrealistic density profile at z ~ 0. In order to measure the positivity 
of AfZ, it is more useful to consider the accumulated A t z [35]. 
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Figure 6. A(Atz) / A^z using sky-averaged analysis for 5 yr observations using DECIGO/BBO 
only (solid), and with joint searches of DECIGO/BBO and ET (dashed). 

In this paper, we have assumed that the correct model is ACDM but it can be other dark 
energy model. However, it has been shown that a variety of dark energy models predict the 
similar values of Atz to ACDM one (see Quartin and Amendola [17] and references therein). 
Therefore, we emphasize that our results are not restricted to ACDM only. 

Unfortunately, it seems very difficult to measure Atz with the ground-based detectors, even 
with the third- generation ones like ET [36]. (The noise curve is shown as dotted curve in Fig. 1.) 
The main reasons for this is because it has less number of GW cycles, the event rate is smaller, 
and the frequency is higher which leads to smaller contribution of "-4PN" effect. However, ET 
can help improving the measurement accuracy of Atz when combined with DECIGO/BBO. 
Figure 6 shows the results comparing the measurement accuracy of Atz with DECIGO/BBO 
observation only (solid), and the one obtained by the joint search of DECIGO/BBO and ET 
(dashed.) We see that the errors are slightly reduced for the latter case. 

There can be peculiar acceleration of each binary which acts as an additional "noise" 
when measuring the cosmic acceleration. However, Amendola et al [37] have estimated the 
peculiar accelerations for typical clusters and galaxies, and found that they are almost the same 
magnitude as the cosmological acceleration. While the peculiar acceleration is different for each 
binary source, the cosmic acceleration is a universal effect. As such, we can safely neglect the 
effect of peculiar accelerations as a noise source since the one from the detector overwhelms it. 
(See also Uzan et al [38] for the discussion about the peculiar acceleration.) 

In this paper, we have assumed that the binary orbits are circular. Once we include the 
eccentricity into parameters, this may have large degeneracy with Atz. We need to estimate 
how the measurement accuracy of Atz is reduced for the eccentric binaries in future. Also, we 
note here that we have used the LTB metric as a simple effective model. To be more realistic, we 
need to use more sophisticated models such as Swiss-Cheese model [39] (see a recent discussion 
by Celerier [40]). We leave these issues for future work. 
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Appendix A. Binary gravitational waveform without the effect of cosmic 
acceleration 

The restricted binary gravitational waveform in Fourier space without cosmic acceleration is 
given as [31, 41] 
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where the amplitude A is given as A = ^^ 2/3 ^jfj— and the phase up to second Post-Newtonian 



(PN) order can be expressed as [21, 31] 
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Here. 



is the coalescence phase and A po \ ia , <£ po i,a anci ¥>d are given as [29,31] 
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where L is the unit vector parallel to the orbital angular momentum, N is the unit vector 
pointing towards the center of mass of the binary, T=lyr and i?=lAU. and F* are the 
beam pattern functions defined as 



F*(0 s ,<h,il>s) 
F+(6 s ,4>s^s) 

F*(e s ,<t>s,i>s) 



-(1 + cos 2 6s) cos(20 s ) cos(2^s) - cos(6» s ) sin(2^s) sin(2^s), (A.7) 

-(1 + cos 2 S ) cos(20 s ) sin(2</> s ) + cos(# s ) sin(2^ s ) cos(2^s), (A.8) 

F+(fl s ,0s-^/4,Vs), (A.9) 

^i x (0s,0s-vr/4,^ s ), (A.10) 



where the direction (6s, 4>s) is measured in the detector's frame and tps is the polarization angle. 
These angles and L ■ N are related to the angles (6s,4>s,Gl,4>l) measured from the center of 
mass of the solar system as explained in Appendix A of Ref . [29] . We have neglected the spin- 
spin coupling at 2PN order since it has been shown that this effect is negligible for NS/NS 
binaries [32]. t(f) in Eq. (A.l) up to 2PN order is given as 
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with t N (f) = (5/256)M z (irM z f)- 8 / 3 2 . When we take the average over the direction and 
orientation of sources, Eq. (A.l) just reduces to 



^(/)lio k a;c7 raged) = ^W- 7 /V*(/). (A.12) 
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2 When we consider the effect of cosmic acceleration, there appears a leading "-4PN" correction in t(f), but we 
have checked that it does not affect our results. 



